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ABSTRACT 

We present the first fully calibrated H 2 1-0 S(l) image of the entire 30 Doradus nebula. The 
observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 
4-meter Blanco Telescope. Together with a NEWFIRM Bry image of 30 Doradus, our data reveal 
the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H 2 - 
emitting area, which extends from the northeast to the southwest of R136, is a photodissociation 
region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus 
are photodissociation regions viewed edge-on. Based on the morphologies of H 2 , Bry , 12 CO , and 
8 fim emission, the H 2 to Bry line ratio and Cloudy models, we find that the H 2 emission is formed 
inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the H II region, 
in a relatively low-density environment < 10 4 cm' 3 . Comparisons with Bry , 8 /im, and CO emission 
indicate that H 2 emission is due to fluorescence, and provide no evidence for shock excited emission 
of this line. 


1. INTRODUCTION 


Starburst feedback is vital in galaxy evolution, as 
it is important in un b inding large mo l ecular clouds 
(iKrumholz et al.l 120061 lFa.ll et aLl 12010b Murray et al j 
I2010h . driving gravitational collapse inside molecu¬ 
lar clouds and triggering sequential star formation 
(lElmeereen fe Ladall977UOev et aLll2005l:IZavagno et al.l 

I2010H. driving turbulent motions within the clouds 
(|Matzneil 120021! . eroding molecular clouds by photo¬ 
evaporat i on (Whitworthlll979l iWilliams fc McKeall997l : 


iMatzneil I2002H . and determining emission line spec¬ 
tra of the photoionized regions dBinette et al.1 119971: 
Dopita et all 120051 120061: iDrainel 12011: lYeh fe Ma.tzueil 


2012t lYeh et al.l 120131 : iVerdolini et, al.ll2013fl . While the 


effects of massive star feedback have been extensively dis¬ 
cussed in the literature, a critical piece of information is 
missing: how does energy and momentum feedback from 
massive stars affect the molecular clouds’ physical prop¬ 
erties? We focus on the spatial distribution of molecular 
and ionized hydrogen emission, signals which highlight 
ionization fronts and working surfaces in this region of 
active stellar feedback. 

A prime site for the origin of H 2 emission is at the 
boundary of ionized gas, shocks, and molecular clouds, 
which is the working surface for many forms of stellar 
feedback into the dense gas. Because the excitation 
of H 2 molecules is sensitive to the density structure, 
radiation hardness, or mechanical energy input from 
shocks, H 2 ro-vibrational transitions are unique tracers 
to probe physical properties of the ISM under massive 
star feedback. The 30 Doradus nebula (30 Dor) in the 
LMC is a well-studied and one of the nearest (50 kpc) 
starburst regions, at which distance a very high spatial 
linear resolution (1" = 0.2 pc) can be achieved. The 
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30 Dor nebula is dominated by the very young star 
cluster R136, which produces 1 0 516 hy droge n-ionizing 
photons per second dCrowther fe Dessartlll998P . ionizing 
neutral material an d driving it outward . Although the 
Dragonfish Nebula ( Rahman et al1l201lD and NGC 3603 
dConti & Crowtherl 120041 ! are much closer (< 10 kpc) 
and produce hydrogen-ionizng luminosities of 10 51 ' 8 and 
10 515 s _1 , respectively, 30 Dor’s convenient location 
out of the galactic plane allows a straightforward 
comparison to more distant starburst regions. Many 
observations of 30 Dor have been presen ted i n the wave- 


lengths from the 3 

brav to radio 

Townslev et al. 

2006 

; Clni & Kcnnicutt 1994 

Poglitsch et al. 

1995 

: Rubio et al.l 

1998g IndebetoLi 

w et al.l 200(1 

Johansson et al. 199S 

Indebetouw et al.l 

20lll), however, 


a fully calibrated H 2 map of the entire 30 Dor has never 
been produced. Wide field of view images of multiple 
ISM components are critical for followup studies on 
dynamics of the region, and for the selection of sites 
for targeted observations, such as high-resolution spec¬ 
troscopy and interferometric studies, which necessarily 
involves smaller fields of view. 

Observations of H 2 in parts of 30 Dor have been done 
by iPoglitsch et al.l (|1995f) (P95) and iRubio et al.1 (|1998H 
(R98). P95 reported that the H 2 1-0 S(l) morphology 
appeared fragmented (~1", 0.2 pc clumps), and sug¬ 
gested that the emission originated from dense molec¬ 
ular clumps. The P95 data is observed in a 3'x3' area 
in 30 Dor, which is a very small fraction of the region, 
while the R98 H 2 data were collected in a small area (< 
2'x2') and not calibrated. These observations were lim¬ 
ited in sensitivity and field of view by the instruments 
available at the time, and the photodissociation region 
(PDR) physical conditions derived from P95 and R98 do 
not fully represent that in the entire 30 Dor. 

In this paper, we present the first fully calibrated 
H 2 1-0 S(l) image as well as a Bry image of the full 
nebula, with a 1.0" angular resolution. We describe 
the observations and data reduction in § ei in § m we 
show the H 2 and Bry morphologies and determine line 
ratios, identify areas of interest for further analysis, and 
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investigate the spatial relations between H 2 and Br 7 . 
In § [4j we present photoionization models using Cloudy, 
in order to constrain the range of physical parameters 
inside the 30 Dor PDRs by comparing the modeled and 
observed H 2 to Bry line ratio, and to explore the issue 
of bright line contamination we discovered during data 
reduction in the H 2 image. The origin of the H 2 1 0 
S(l) emission is discussed in § [5] Finally, we summarize 
the paper in § [6] In the paper, all H 2 emission refers 
to the H 2 1-0 S(l) transition, unless indicated otherwise. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed 30 Dor using the NO AO Extremely Wide - 
Field Infrared Imager fNEWFIRM: iProbst et all (|2008T ) 1 
on the CTIO 4-meter Blanco Telescope, over three half¬ 
nights on November 10, 11, and 12, 2010. NEWFIRM 
has a field of view of 28 , x28' 1 and its pixel scale is 0.4" 
per pixel. The H 2 1-0 S(l) (2.12 pm) and Bry (2.17 pm) 
emission line data were taken using the 2124 nm H 2 and 
2168 nm Bry narrow band filters, respectively. See Ta¬ 
ble Q] for filter parameters. The broad-band continuum 
data were collected using the Ks filter. The total expo¬ 
sure time of the H 2 image was 210 minutes, and 14 min¬ 
utes for the Bry image. The observations were dithered 
in a random pattern in a 30" box to fill the gaps between 
detector arrays. The photometric standard star S121-E 
was observed in both 2124 nm and 2168 nm filters to 
serve as the flux calibrator. Table [I] summarizes details 
of the observations. Because the angular size of 30 Dor 
in the sky is about the same as the size of NEWFIRM 
field of view, in order to obtain sky images free of neb¬ 
ular emissions, we nodded the telescope on and ofl0 the 
target following the sky-target-target-sky sequence. 

Data reducti on was carried out u sing the NEWFIRM 
pipeline VI.3 (jSwaters et all [20091 ). Dark subtraction 
was first applied to the data, followed by a linearity 
correction and flat fielding. Sky background was 
determined by taking the median of four preceding 
and four subsequent off-source sky exposures which 
are free of extended emission, and the background 
level was scaled to match that in target images and 
then subtracted. An astrometric solution was obtained 
using the 2MASS catalogue, and all images were re¬ 
projected and stacked. The sky background was then 
redetermined and subtracted by masking objects (in¬ 
cluding stars and extended nebular emission) detected 
in the first pass, and new stacks of images were produced. 


2.1. PSF Matching 

We carried out a photometric analy sis of field stars us¬ 
ing t he software package SExtractor dBertin fe Arnouti 
fl996lf . Stars with detection higher than 5 a in areas 
free of extended nebular emission are identified by 
SExtractor, and their photometric parameters, such as 
flux and FWHM, are recorded in catalogs. We first 
extracted the mode of the distribution of FWHM in each 
stacked image, namely Bry , H 2 , and Ks images, as the 
representative seeing in each image. It showed that Bry 
and Ks images have better seeing (< 1.0") than the H 2 
image. We then convolved the Bry and Ks images with 

4 The telescope was nodded off-target 1 degree north of R136. 


gaussian kernels until the convolved FWHM matches 
that of the H 2 image. After the gaussian convolution, 
the seeing in H 2 , Bry , and Ks images is 1.0", and this 
is the resolution in every image presented in the paper 
unless indicated otherwise. 


2.2. Flux Calibration 

The standard star S121-E dPersson et al.lll998tl was ob¬ 
served as a flux standard. The Bry and H 2 filters are 
centered closely at the Ks filter central wavelength, and 
we apply the magnitude-to-flux density conversion fac¬ 
tors derived from the S121-E data, to flux calibrate the 
30 Dor data. Data were taken under stable weather and 
nearly constant airmass, therefore the major uncertain¬ 
ties of measured fluxes come from bright emission line 
contamination in the Ks filter, which are discussed in 

§mi 

2.3. Continuum Subtraction 

To produce emission line images, one must subtract 
continuum emission in the narrowband (NB) data. The 
ideal way to subtract continuum in a NB image is to use 
an off-line center NB filter with identical FWHM, and 
such NB image shows continuum without contamination 
of any emission lines. However, off-line NB filters were 
not available when the observations were carried out, we 
therefore carried out continuum subtraction using the 
broadband (BB) Ks data. 

Mathematically, the observed fluxes in each filter 
(2.12 /rm, 2.16 pm, and Ks) can be expressed as 

T 2124 Fb 2 V FA 2124 , 

-K2168 — Fsr'y V F c t 21G8 , 

Fks — F h , 2 T FB r ry T F u k T F c t Ks , 

where F 2 i 2 4 , F 2168 J and Fx s are fluxes measured in 
the 2.12 pm, 2.16 pm, and Ks filters, respectively. The 
convention ‘ct’ labels the continuum emission fluxes mea¬ 
sured in a filter, and ! uk’ marks unknown emission line 
fluxes contained within the Ks filter. 

Continuum subtraction in Bry and H 2 images then 
will produce: 

Tiugg rxFj^ s — Fj^ r ry -f- F ct2les 

cr(F// 2 T FB r ry T F u k T F c i Kb ) , 

yielding 

„ -F 2168 — a{Fxs — Fh 2 — F u k ) . . 

r Bri — Z ) (1J 

1 — a 

and 


-F2124 — @ Fk s 

yielding 

Fh 2 


h'ff'2 V h’cf 2l 24 

P(Fh 2 + Fbt-1 + Fuk + F ctKs ) , 
T2124 — P(F Ks — Fb^) 

. W 




where aF Ks = F ct2168 and f3F Ks = F ct2124 . The equa¬ 
tions are applicable at every pixel in the image. The 
values of a and /? in principle should be close to the 
FWHM ratio of the Bry to Ks filter and H 2 to Ks filter, 
respectively. 
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Table 1 

30 Doradus Observation Summary 


Line 

Vacuum Wavelength 

Filter FWHM (nm) 

Filter name 

Observed Dates a 

Total Integration Time (minutes) 

Ks continuum 


320.0 

Ks 

Nov 10, 12 

22 

H 2 1-0 S(l) 

2.121 fim 

24.0 

2124 nm 

Nov 10, 11, 12 

210 

Bry 

2.166 fim 

24.4 

2168 nm 

Nov 11 

14 


a Data were taken in year 2010. 


We determine a and /? empirically by evaluating the 
stellar flux ratio of 2.16 /jm to Ks, and 2.12 /mi to Ks. 
We employed SExtractor to extract stellar fluxes in the 
2.12 /im, 2.16 ;im, and Ks images. The scaling factor a 
and /3 are 0.07 and 0.08, respectively, which is in good 
agreement with the filter FWHM ratios. The above equa¬ 
tions are correct when the CCD response is linear in all 
three filters. However, stars with counts > 10,000 Ana¬ 
log/Digital Units are saturated, i.e. the CCD response 
becomes nonlinear, and they cannot be completely sub¬ 
tracted. Therefore we exclude these stars in the analysis. 

We found that the BB Ks data contain Bry , H 2 , and 
possibly other emission lines, as well as the continuum 
emission, which introduce contamination. We evaluate 
and discuss the bright line contamination issue in § 12.3.11 


2.3.1. Bright Emission Line Contamination 

No direct information is available to us on the possi¬ 
ble emission lines that, other than Bry and H 2 , might 
contribute to the Ks filter emission. However, He I is a 
likely candidate since it is observed in other star forming 
regions. For example, heliu m emission lines are rep orted 
in M16 in the 2 /tm regime (jLevenson et, al.l [200(1 here¬ 
after LOO), in addition to Bry and H 2 lines. Amongst 
the detected emission lines, Bry is the brightest in M16, 
and a bright He I line at 2.06 /rm is 70% of the total flux 
of Bry . If the He I -to- Bry ratio is the same in 30 Dor 
as in M16, the contamination from He I in the Ks filter 
will be noticeable. 

Bry emission line in 30 Dor is likely the brightest 
amongst the emission lines in the Ks filter. If the He I 
line has the same relative strength in 30 Dor as in M16 
(70%), Equation[l]shows that it is equivalent to about 5% 
of the continuum-subtracted Bry emission, assuming it 
is distributed in the same way. Therefore continuum sub¬ 
traction in the Bry image may not be severely affected 
by bright emission line contamination, and Equation |T) 
can be approximated as 


F'bt'y — 


E2168 — uFks 
1 — a 


( 3 ) 


In the case of H 2 continuum subtraction, however, 
bright emission line contamination becomes significant. 
Following Equation [2j the first-pass continuum sub¬ 
tracted H 2 shows that if no correction was made, then 
there would be a noticeable over-subtracted area. In¬ 
deed, without correction, we noticed a region of negative 
emission representing the Bry emission in the vicinity of 
the nebula. With careful visual inspection, we found that 
the negative components well resemble the morphology 
of the highest surface brightness Bry emission. We then 
inspected the first-pass continuum subtracted H 2 and 
Bry images on the pixel-to-pixel basis, comparing the 
pixel values of the negative component in the H 2 image 


to the Bry image. There is a tight correlation between 
the negative H 2 pixels and brightest Bry pixels, confirm¬ 
ing that the majority of over-subtraction comes from line 
emission which correlates strongly with Bry emission. 
Let 


T P , _ T2124 — P ( Fks 

= F u k , 

Equation [2] can be rearranged as 


F'bt'y) 


and 


F’h 2 =Fh 2 


P 

1-/3 


^FBr 7 


( 4 ) 


The slope in Equation Q] is -0.10, evaluated from 

fitting the residual data. Therefore we correct continuum 
subtraction of the H 2 image by 

F' H2 =Fh 2 — (—0.10 )F Bri 

We inspected the pixel values of the corrected H 2 im¬ 
age F' h , 2 , and the majority of the negative pixels after the 
correction have values around 0. In fact the factor 0.10 
is higher than the expected value /3/(l-/3)= 0.087, which 
indicates that we have not only corrected for the con¬ 
tamination introduced by the Bry emission, but also the 
contamination from other unknown bright emission lines. 
The bright-line contamination is thus (0.10-0.087)/0.087 
= 0.15, or 15% of the Bry emission flux, which corre¬ 
sponds to 1.3% of the continuum flux. 

We found that, after correcting the H 2 image for the 
contamination, some negative pixels still persist in the 
areas very close to R136. Those pixels can only be cor¬ 
rected by = 0.14, instead of 0.10 in Equation [I] 

However this leads to over-correction in the image, i.e. 
bright Bry features become prominent in the H 2 im¬ 
age, which indicates that the contaminating emission 
is stronger relative to Bry near the cluster than fur¬ 
ther away. The corresponding total brightness of this 
contamination is 61% of the Bry flux, assuming that 
the emission line flux is distributed in the Ks band in 
the same way as Bry , which leads to additional 4% of 
the continuum flux, which is also the uncertainty in the 
contamination-corrected H 2 image. We suspect that the 
He I line at 2.06 /im is the major source of continuum 
contamination other than Bry , and its contamination 
becomes more significant in the central region of 30 Dor. 
Several He I lines are detected in the M16 Hll region, 
in addition to Bry and H 2 1-0 S(l) (LOO). As noted 
earlier, the brightest He emission line detected in M16 
in the K band is He I at 2.06 /im, and its flux of the He I 
line is 70% of the total flux of Bry . The 61% of the 
Bry flux contamination we have empirically estimated 
for the region near R136 is as significant as that in the 
M16 case. Note that variations in the continuum slope, 
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such as those induced by variations in extinction, will 
also not be consistent with constant values of our a and 
f3 parameters. 

We do not have any He line data in 30 Doradus to 
constrain the degree of contamination, therefore we turn 
to Cloudy simulations to explore this issue, which is dis¬ 
cussed in §2) In this paper, the H 2 image is empirically 
corrected for contamination following Equation [4] 

3. RESULTS 

3.1. H 2 and Br^ Morphologies 

Fully calibrated H 2 and Bry images of 30 Dor are pre¬ 
sented in a 3-color composite image in Figure [T] Red is 
H 2 , blue is Bry , and green is the Ks band continuum. 
Both Bry and Ks are stretched logarithmically to em¬ 
phasize areas of highest surface brightness, while H 2 is 
displayed in the linear scale because the line brightness 
dynamic range is much smaller. The star cluster R136 is 
marked by a black circle. Note that the pulsar wind neb¬ 
ula N157B, which is located at the same distance and T 
to the southwest of 30 Dor, is also present in the image. 
In this paper we focus only on the analysis of 30 Dor and 
will ignore N157B. 

The Bry emission reveals the spatial distribution of 
ionized gas in 30 Dor. Areas with highest surface bright¬ 
ness appear to trace an arch structure extending from 
northeast to southwest of the R136 cluster, consistent 
with the ionized gas morphology reported in the lit¬ 
erature dChu fc Kennicuttl 11994 iPoglitsch et all 119951 
[Pellegrini et all 120101 1. To the north and west of R136, 
Bry appears quite filamentary and its surface brightness 
becomes lower. To the southeast of 30 Dor, the Bry 
morphology reveals multiple shell structures, enveloping 
pillars and clumpy features. The total flux in an area 
3.1'x 3.7' is measured as I.02xl0 -1 ° ergs -1 cm -2 . Al¬ 
though this is higher than the total Bry flux reported in 
P95 (4xl0 -11 ergs -1 cm -2 ) measured in the same area, 
our Bry image detects fainter structure than that in the 
P95 image. 

The most prominent H 2 emission is seen to the north¬ 
east of R136, in conjunction with the bright Bry arch and 
coincident with lower surface brightness Bry emission. 
This prominen t H 2 emitting area, which is also known as 
30Dor-10_[e.g. iJohansson et all [l998l Hiidcbctou w et al.l 
l20Ml20m also marked as Area A in Fig. [I| spans at 
least 3' by 3' in the plane of the sky, which corresponds 
to 36x36 pc at a distance of 50 kpc. The H 2 emission 
in 30Dor-10 seems somewhat disordered and extended, 
with clumps close to R136 and filaments extending away 
from R136. The total H 2 flux measured in the area is 
5.16xl0 -12 ergs -1 cm -2 after masking bright saturated 
stars. The H 2 emission to the west of R136 appears 
much clumpier and mixed with high surface brightness 
Bry emission. To the north and northeast of R136 (Area 
B), H 2 appears in the form of filaments which seem to 
form a chimney pointing away from R136, and their mor¬ 
phology is poorly correlated with the Bry filaments in 
the same area. Prominent pillar features are seen in the 
southeast of R136, pointing towards the ionizing source 
R136 and are encompassed by the Bry emission, such as 
that in Area C. These Bry envelopes have sharp outer 
boundaries, with radii of curvature significantly greater 
than those of the H 2 pillars they envelop. This is sugges- 



Figure 1. Three-color composite image of 30 Doradus. North is 
up and East is to the left. Red: H 2 , displayed in linear scale; green: 
Ks continuum in logarithm stretch; blue: Br 7 in logarithm stretch. 
R136 is marked by a black circle. The seeing is ] .0". H 2 -emitting 
areas of interest are marked by white dashed ellipses, including 
(A) the bright H 2 band to the northeast of R136, (B) filaments 
pointing north-ward of R136, and (C) pillars to the southeast of 
R136. A zoomed-in figure of Area C is shown on the right. A 
bright H 2 finger in Area A is indicated by a black arrow. 

five of a photo-evaporative flow bounded by the pressure 
of hot gas. No H 2 emission was detected to the north¬ 
west of R136. Overall, the Bry and H 2 emission appear 
to form walls of the cavities or holes seen in 30 Dor. The 
observed H 2 emission is located well within the nebula 
(defined by optical broadband dat a). We will discuss 
Area A, B, and C in detail in § 13.21 

3.2. Areas of Interest 

We identified areas of interest for further analysis by 
comparing the H 2 , Bry , and 12 CO morphologies as 
shown in Figure [2] and Figure |T] The H 2 and Bry 
data were superim posed with the 12 CO 1-0 data of 30 
Dor (iWong et al.l 1201 ljl 5 ! . The Bry image in Figure [2] 
is displayed in grayscale to guide the reader’s eyes on 
the ionized gas morphology, while the H 2 contours are 
convolved to a 4.5" resolution. We found that bright H 2 - 
emitting areas spatially correlate very well with the bulk 
of CO emission (Figure [2j). Based on the morphological 
correlation between H 2 and CO emission, and H 2 and 
Bry emission, three areas are identified as areas of inter¬ 
est: (1) Area A: the northeastern band of H 2 emission 
with highest H 2 surface brightness, and it spatially coin¬ 
cides with the brightest CO emission in 30 Dor. This area 
is also known as 30Dor-10 as noted earlier; (2) Area B: 
the filament pointing north and away from R136, which 
is relatively bright and spatially correlated well with the 
peak of CO emission in the same area; and (3) Area C: 
pillars to the southeast of R136 at the outer shell of 30 
Dor, which have clearly defined morphology in the H 2 
emission and are surrounded by Bry-emitting envelopes; 
no CO emission is detected in this area (by the ATNF 

5 We cha se f he IWong et. al. 120111) CO data for comparison in¬ 
stead of the llndebetouw et aH (1201 3b ALMA data, because we are 
interested in the CO distribution in the entire 30 Dor. The ALMA 
image field of view is too small, and the interferometric data is not 
sensitive to extended structures. Nevertheless, the ALMA data 
will be useful for detailed studies of smaller objects of interest in 
the region. 
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Mopra Telescope. IWong et al.1 (|2011[ V). A zoomed-in fig¬ 
ure of Area C is shown in Figure [l] 

H 2 was detected in an earlier observation (P95) in 
the western part of Area A, in a much smaller 2'x 2' 
area. The H 2 morphology observed by P95 was reported 
very clumpy, which is consistent with the H 2 morphology 
shown here (Figure [1]). The peak H 2 surface brightness 
in the P95 result is 5.41 xlO -15 erg s -1 cm -2 arcsec -2 . 
However, this peak intensity arises from a bright star 
in the field thus may not be reliable. In our wide- 
held observations, more H 2 emission is seen towards 
the eastern part of Area A. After excluding bright satu¬ 
rated stars, the maximum H 2 surface brightness is mea¬ 
sured 2.15xl0 -15 ergs -1 cnr -2 arcsec -2 . H 2 in Area A 
al so co incides with the bulk of CO emission in 30 Dor 
(IWong et al J 1201 ll) (Figure [2]). No H 2 emission is seen 
completely uncorrelated with either the Bry or 12 CO 2- 
1 emission. With only one ro-vibrational transition of 
H 2 emission, and without making further assumptions, 
we estimated mass contained in the upper state of the 
observed H 2 transition in Area A to be 0.01 Mg, esti¬ 
mated with the Einstein A coefficient Ag( 1 )=2.09x 10 - ' 
s -1 (|Turner et alJll977t IWolniewicz et al.l 11998ft . 0 This 
mass is orders of magnitude lower than the molecu- 
lar mass estimated from CO emissions of 8.5xlO 4 M Q 
(|Pineda et al.l 120121 : I Johansson et al.l Il998l) . which im¬ 
plies that the observed H 2 emission requires the exis¬ 
tence of only a negligible fraction of the total molecular 
mass. 

Area B also shows high H 2 intensity, and the 
maximum H 2 surface brightness in this area is 
6.44xlO -16 ergs -1 cm -2 arcsec -2 . This region was out¬ 
side the field of view of the P95 observations. Area B 
lines up well with its CO counterpart (Figure [2]), similar 
to Area A. 

Area C was outside the field of view of the P95 ob¬ 
servations as well. The peak surface brightness is mea¬ 
sured 1.76x 10 -15 ergs -1 cm -2 arcsec -2 , comparable to 
that measured in other areas. Although Area C does 
not appear to have any CO counterpart, its H 2 emis¬ 
sion indicates the presence of molecular clouds in that 
area. The morphology of the pillars is very similar to 
t he ones observed in smaller Hll regions, such as M16 
dHester et allI1996L LOO). We highlight few interesting 
morphological features in Area C, which will lead to fol¬ 
lowup detailed studies of the region: (1) the H 2 emis¬ 
sion is enclosed within the Bry emission which extends 
toward R136, (2) the H 2 emission appears at the leading 
edge of regions where other molecular cloud tracers are 
present, such as the 8 //m emission, and (3) the Bry emis¬ 
sion is bounded by a sharp outer boundary whose radius 
of curvature is significantly larger than that of the H 2 , 
which implies possible pressure confinement by hot gas. 
A crude estimate of the projected distance between the 
Bry envelopes and the pillars is ~ 3 pc. Area C’s neigh¬ 
bor clumps also show similar separation between Bry 
and H 2 . We will return to discuss this point in § 15.11 

3.3. The Spatial Relationships Between Ionized and 
Molecular Gas in 30 Doradus 

6 Note that this mass estimate yields a lower limit of molecular 
mass, because it does not include the mass associated with all other 
levels. 



Figure 2. CO 1-0 contours (white contours, IWong et al.l!201ll 'l 
superimposed on the 3-color composite image as shown in Figure 
1. The CO contours start at 10% of the peak intensity (12.4 K km 
s — and increase in 10% steps. The CO data angular resolution 
is 60". R136 is marked by a black circle. 

The overall morphology of the observed H 2 emission is 
well correlated with the Bry distribution within 30 Dor - 
when Bry emission is present, H 2 emission is either seen 
very nearby the Bry emission, such as the outer shells 
and pillar features, or appears to be a counterpart of the 
Bry emission, such as Area A and the region just west of 
R136. 

Bright H 2 -emission is seen to be spatially coincident 
with the Bry emission in Area A (Figure [I]), except for 
an excess of Bry emission around the limb of the “fin¬ 
ger” seen in H 2 . This suggests that we are viewing the 
H II region roughly face-on, except around the limb. Fur¬ 
thermore, the CO-emitting gas in the core of the finger 
is likely to be marginally optically thick at 2.2 /mi, judg¬ 
ing from the ALMA obs ervations of 13 CO presented by 
llndebetouw et al.l (120131 ! . They estimate 13 CO column 
densities > 5 x 10 15 cm -2 along this finger (and clumps 
an order of magnitude higher column; their Fig. 10), and 
adopt Ah 2 = 5 x 10 6 Ah3co; along with a dust opac¬ 
ity of roughly 1.3 x 1 0 -2O cm 2 per H atom at 2.17/tm 
l)Weingartner fe Drainel[200il l. we estimate the charac¬ 
teristic optical depth of the finger to be ~ 0.65. It is 
therefore quite likely that most of the observed emission 
is emitted on the Earthward side of the molecular finger. 
The absence of obvious shadows associated with high- 
column clumps of 13 CO corroborates this interpretation. 

On a similar basis, we infer that emission seen just 
west of R136 is likely to originate from behind the cl uster 
as se en from Earth), on the basis that the IWong et all 
201 If ! CO contours extend into this region. This ge¬ 
ometry is furt her supported by recen t optical emission 
line studies bv iPellegrini et ai] (12011IK P111. With opti¬ 
cal emission line ratios together with modeled ionization 
parameters, Pll suggests that Area A and the region 
just west of 30 Dor are located 60 pc behind R136. 

Area B may share the same geometry as the regions we 
discussed above - this filament of molecular gas should 
be behind the Bry emission and R136. Although the 
Bry morphology is filamentary with rather poor corre¬ 
lation with the H 2 emission, Bry is seen to partially 
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coincide with the H 2 and CO emission in the area, in¬ 
stead of being on the edge of molecular clouds. This sug¬ 
gests that the Br 7 emission is not shielded by molecular 
clouds, and we are viewing this region face-on. Lacking 
high-resolution CO data, however, we cannot draw this 
conclusion with certainty. 

Area C is likely oriented differently relative to 30 Dor 
than are Areas A and B. The bright Bry emission in this 
region shows relatively poor spatial correlation with the 
H 2 emission. Bry is seen closer to R136 in the projected 
plane of sky, and it envelops the H 2 emission, which has 
a relatively small projected area. Although Bry extends 
outward from the H 2 emission (and toward 30 Dor) to a 
similar degree here as in Areas A and B, the much smaller 
H 2 -emitting regions give it more of the appearance of an 
evaporation flow. 

In all three regions there are locations where Bry emis¬ 
sion extends away from a boundary of the H 2 emission, 
generally in the direction toward 30 Dor, before terminat¬ 
ing at an enclosing boundary about two to three parsecs 
away. To be specific, in Figure [T| we measure projected 
standoff distances of 2.5 pc perpendicular to the finger in 
Area A (to the southeast); of 2.2 pc (toward 30 Dor) from 
the small H 2 blob at the southern end of Area A; of about 
1.9 pc (toward 30 Dor) from parts of the H 2 ridge in Area 
B; and of 1.9 pc (toward 30 Dor) and 2.2 pc (perpendic¬ 
ularly to the northeast) from the H 2 pillars in Area C. 
All of these characteristics - the tendency of Bry emis¬ 
sion to enclose H 2 , its tendency to extend a few parsecs 
toward 30 Dor, and its tendency to meet a sharp bound¬ 
ary, are expected to reflect the physical origin of these 
two types of emission, as well as the relative importance 
of effects such as photo-evaporation and confinement by 
stellar wind pressure. 

3.4. H 2 to Bry Line Ratio 

The morphological relations between H 2 , Bry , and 
CO emission are indicative that the observed H 2 emis¬ 
sion traces the PDRs in 30 Dor. With data of only one 
H 2 emission line, we do not have sufficient information 
to firmly constrain the physical quantities in the PDRs. 
Nevertheless, the line ratio of H 2 to Bry is a useful 
guide to delineate the spatial distribution and structure 
of molecular gas relative to ionized gas in the PDRs, es¬ 
pecially when a PDR is viewed face-on. Combining with 
numerical modeling efforts, the observed H 2 to Bry line 
ratio offers a hint of the physical properties of molecular 
and ionized gas in 30 Dor. 

A line ratio map H 2 to Bry is shown in Figure[3] Pixels 
in the Bry image were clipped at a 3<r level, and pixels 
with S/N higher than 50qJ detection in the H 2 image 
were masked in order to exclude bright saturated stars. 
A line ratio map then was convolved to a 4.5" resolution 
with a gaussian kernel. 

Overall, areas with higher H 2 / Bry ratios in 30 Dor 
are clumpy with ratios of 0.2 to 0.5. Most areas show 
lower H 2 / Bry ratio ~0.1, which agrees with the same 
line ratio observed in M16 (LOO, derived from the total 
fluxes). Higher H 2 / Bry ratios are seen in Area A, Area 
B, filaments north of Area B, and some isolated pillar 
features (including Area C) at the outer shells of 30 Dor. 

7 An empirical value to mask saturated bright stars. 



Figure 3. H 2 to Br 7 line ratio map of 30 Dor, convolved by 
Gaussian kernels to a 4.5" resolution. The colorbar marks the H 2 
to Bry ratio; R136 is marked by a cyan cross. 

In Area A, the H 2 to Bry ratio across the area ap¬ 
pears clumpy, with some localized high line ratio areas 
and ‘voids’. The maximum ratio is 0.5, and the overall 
ratio is > 0.2. The clumpy H 2 distribution and H 2 to 
Bry line ratios seen in Area A indicate that FUV radi¬ 
ation could penetrate deeper into the molecular clouds. 
The maximum line ratio at Area B is 0.45, and the over¬ 
all radio is ~ 0.3. The line ratio morphology suggests 
that Bry across this area is fairly filamentary. Area C 
displayed a high H 2 to Bry line ratio of 0.3, which coin¬ 
cides with the bright ^-emitting area at the tip of the 
pillars. 

We note that the two ISM components traced by Bry 
and H 2 arise in adjacent but noticeably separate areas, 
and the local, apparent line ratios are subject to signif¬ 
icant projection effects. The projection effect is most 
severe in regions where the ionized gas envelops a molec¬ 
ular pillar viewed edge-on, such as in Area C. 

4. PHOTOIONIZATION MODELS 
4.1. Motivations 

Our main motivation to carry out numerical simula¬ 
tions is to explore a wide range of physical conditions in 
H II regions and PDRs, and obtain emission line intensi¬ 
ties as a function of physical conditions. The observed 
Bry and one single H 2 emission provide very limited in¬ 
formation, but with numerical simulations, we are still 
able to learn something useful about the physical con¬ 
ditions in 30 Dor. Another motivation is to explore the 
issue of bright emission line contamination in the data. 
A series of H 2 emission lines, hydrogen recombination 
lines, and He recombination lines are present in the 2 pm 
regime. We found evidence of additional emission line 
contamination (other than Bry ) in the Ks continuum, 
affecting continuum subtraction in the H 2 image (see 
§ 12.3.11) . We have empirically corrected for the contami¬ 
nation, and the estimated contamination level is in good 
agreement with similar type of observations in other H II 
regions. Numerical simulations are helpful to investigate 
feasible origins of contamination and consequences. 

4.2. Model Parameters 

We generated models of simple Hll regions to study 
emission line intensities and physical conditions of molec¬ 
ular clouds in 30 Dor. First we used Starburst99 
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Figure 4. H 2 1-0 S(l) to Br 7 line ratios as a function of molec¬ 
ular gas density (n mo i) and distance to the ionizing source. The 
purple line shows H 2 to Br 7 line ratio computed at a distance of 
30 pc, and the line ratio computed at 60 pc is shown in the gray 
line. Observational constrains on H 2 to Br 7 ratios and n mo i in 
Area A and C are marked by blue dotted and cyan dashed lines , 
respectively. 

(iLeitherer et al.1 11999ft to generate ionizing continuum 
spectra of a massive coeval star cluster at 2 Myr age, 
because the age of R136 is < 2 Myr (|de Koter et al.1 
119981: iMassev fc Hunteil 1 1998ft . The star cluster is as¬ 
sumed to be massive enough to fully sample the initial 
mass function, which has exponents -1.3 and -2.3 between 
stellar mass boundaries of 0.1, 0.5, and 120 M 0 . We 
employed the Geneva high mass-loss evolutionary tracks 
with 0.4 solar metallicity. The Geneva high mass-loss 
tracks are optimized for modeling atmosp heres of high 
mass stars and are recommended bv iMaeder fe Mevnetl 
(11994ft . We adopt Pauldrach/Hillier atmospheres and 
the LMC UV line library. The at mosp he res include 
non-LTE a nd line-blanketing effec ts (jSmith et al.l 12002ft 
for O stars (iPauldrach et al.ll2001ft and Wolf-Rayet stars 
(|Hillier fc Milled 1 1998ft . 

Starburst99 output continuum spectra are fed into 
Cloudy O8.OC0 as th e ionizing continuum of each sim¬ 
ulated H II region. IPellegrini et al.1 (120101 12011ft sug¬ 
gested that the inner 15 pc of 30 Dor lacks ionized gas 
and thus molecular c louds. The ion i zation front near 
Area A identified in IPellegrini et al.1 (12010ft is located 
at a projected distance of 30 pc, and materials behind 
the star cluster R136 are at a characteristic distance of 
60 pc. We therefore set the inner radius of the simu¬ 
lated Hll regions at 30 pc and 60 pc with respect to 
the ionizing sourccfl A wide density grid log(nH)=l 
to 5 is incorporated in the calculations. We adopt the 
dust grain size distributions of the LMC ([Pellegrin i et al.1 
1201 It iWeingartner fe Draindl2001ft . and Cloudy’s default 
ISM abundances at 0.4 solar metallicity. Each calcula¬ 
tion stops when the cloud temperature drops to 30 K, 
well inside a PDR. Integrated line luminosities, Bry at 
2.17 /im, He I at 2.06 /im, H 2 1-0 S(l) at 2.12 /im, H 2 
1-0 S(0) at 2.22 ^m, and H 2 2-1 S(l) at 2.25 /rm, are 
calculated and recorded. 

8 Calculat ions were performed with version 08.00 of Cloudy, last 
described bv IFerland et al.l 119981 1. 

9 Regions modeled at a projected dis tance < 10 p c wo uld require 
ionizing sources other than R136 IPellegrini et alH201HL 


5. MODEL RESULTS AND DISCUSSION 

With the model parameters indicated in § 14.21 we first 
modeled H II regions under perfect force balance between 
thermal pressure of the cloud and incident radiation and 
wind pressure. At the distance 30 pc and 60 pc away 
from the ionizing source, ionized gas density in the force- 
balanced H II regions appears constant, however, in the 
PDR region the cloud pressure can become unrealisti¬ 
cally high when the ionized gas density is high. This 
is because our simple models do not include turbulence 
and magnetic field pressure terms, which are important 
in supporting molecular clouds.The modeled regions thus 
form artificially high thermal pressure beyond the ion¬ 
ization front (IF) in order to meet the force balance cri- 
terium. To resolve this problem, we turned to the con¬ 
stant density model. The Bry emission intensities calcu¬ 
lated in both force balance and constant density models 
are consistent to within 1%, therefore it is reassuring 
that the constant density models adequately represent 
the emission line spectra in the ionized gas, which is im¬ 
portant in affecting the FUV spectra entered into PDR. 
We generated tables of emission line intensities as a func¬ 
tion of the density grid and interpolated the observed 
line ratios using such tables, to obtain plausible range of 
physical conditions in 30 Dor. 

5.1. Molecular Gas Density 

We also modeled simple H II regions with parameters 
described in §[4]and computed the Bry and H 2 1-0 S(l) 
emission intensities as a function of gas density logio nn 
= 1 to 5 (FigureSJ). The computed H 2 to Bry line ratios 
at 30 pc and 60 pc are shown in purple and gray lines, 
and the dotted and dashed lines mark observational con¬ 
straints of molecular gas densities and H 2 / Bry ratios in 
Area A and C. The lower limit of molecular gas densities 
is set by the ionized gas densities evaluat e d from the [SII] 
emission line doublet in iPellegrini et al.l (|2010ft near the 
IFs at Area A and C. The ionized gas densities at those 
locations are found to be 10 2 ' 7 cm -3 and 10 1 ' 9 cm -3 , re¬ 
spectively. The upper limit of H 2 to Bry line ratio in 
each area is set by our observations. 

The H 2 to Bry line ratio in most parts of Area A is 
<0.5, implying that the molecular gas density is < 10 4 
cm~ 3 at 30 pc and < 10 3 cm -3 at 60 pc. The maximum 
H 2 to Bry ratio in Area A is ~0.5, and the corresponding 
molecular gas density would be ~10 43 and 10 31 cm' 3 
at a distance 30 pc and 60pc, respectively. The simple 
comparison of observed and modeled line ratios suggests 
that H 2 in Area A is formed in relatively low density 
areas in the PDR of 30 Dor, close to the surface of the 
molecular cloud in conjunction to the IF, where ionized 
gas density is ~ 10 27 cm -3 . This is consistent with our 
findings in t|3] that the observed H 2 emission arises from 
the PDR near IFs. 

In the Area C pillars, the maximum H 2 to Bry line 
ratio measured is 0.3, where H 2 emission intensity is also 
highest. The projected distance of Area C is greater than 
30 pc, we therefore refer to the 60 pc model. At the upper 
limit of H 2 to Bry ratio of 0.3, the molecular gas density 
would be ^10 2 4 cm -3 . The Area C H 2 emission is seen 
with a projected separation of 3 pc from its Bry enve¬ 
lope. With the ionized gas density 10 1 ' 9 cm -3 , a depth of 
3 pc would corresponds to Ay <0.5. One expects molec- 















































ular hydroge n to form ins ide a m olecular cloud at depths 
Ay ~ 0.13 (|van Dishoeck fe BladdH988f) . which further 
supports the notion that the observed H 2 emission in 
Area C is formed inside the PDR but at a rather shallow 
depth. This conclusion holds just as well for the other 
locations where we have identified Bry extending 2-3 pc 
from an edge in the H 2 emission. As noted in > 13.41 the 
apparent H 2 / Br 7 line ratio determined on small scales 
is subject to projection effects because of this physical 
offset. 

We can compare the observed H ? emissio n fluxes 
with PDR models such 
and 


as 


(1989)(SD89) 


Sternberg fc Dalgarno 


Black & van Dishoeck 


1987|)(BvD87), and constrain molecular gas densi¬ 


ties in Area A, B, a nd C. The far UV radiation field y 
in 30 Dor is ~500 Yn (iPineda et a.l.ll2009t lAndersonetjaL 
12014 where yo (|Drainel I1978H relates to the lHabing 
(|1968li field by a factor of 1.71), while SD89 and BvD87 
predicted that the H 2 1-0 S(l) intensity is <2% of the 
total H 2 intensity at such radiation hardness. Adopting 
the average H 2 emission intensities at Area A, B, and 
C as 2% of the total H 2 intensity, the corresponding 
molecular gas densities in those areas are < 10 4 cm -3 , 
in coarse agreement with our Cloudy calculations. To 
firmly constrain molecular cloud densities in the 30 Dor 
PDR, we will need multiple ro-vibrational transitions 
of H 2 emission in followup studies. With one H 2 
transition and Bry , nevertheless, our analysis suggests 
that molecular clouds associated with the observed H 2 
emission have densities < 10 4 cm -3 . 


5.2. Fluorescence or Shock Excitation? Origin of the 
H 2 1-0 S(l) Line Emission 

NIR molecular hydrogen emission lines in HII regions 
can form either via (1) pure fluorescence excitation or 
(2) shock heating, and the best way to distinguish flu¬ 
orescence from shock excitation is to analyze multiple 
transitions of ro-vibrational molecular hydrogen emission 
lines. Although such spectroscopic data do not yet exist, 
the H 2 morphology, compared with CO morphology and 
H 2 -to- Bry ratios, suggests no evidence of shock excita¬ 
tion. 

As reported in § 13.11 and § E3 the H 2 morphology 
generally correlates well with that of CO. H 2 is likely to 
correlate poorly with CO emission in the case of shock 
excitation, which is often seen in star forming regions 
with active protostellar outflows. A good ex ample is the 
Orion A giant molecular cloud (jDavis et al.ll2009 ). 

The H 2 -to- Bry line ratio is another diagnostic to dis¬ 
tinguish fluorescence from shock excitation. In shock- 
dominated re gions, the line ratio is often found greate r 
than unity ( Puxlev et, al.1 l2000t iMedling et all 120151 ): 
while in ma ssive star-forming regions, the H_> -to- Bry ra¬ 
tio is < 0.6 (l.loseoh et aliri984 iMoorwood fc Oli vail 19881 
iRodrfguez-Ardila et al.l 12004 120051 iRiffel et al.1 1201 011 
The reason behind the distinctive line ratios is fairly 
simple. Shocks cannot excite Bry emission, therefore 
in shock-dominated regions the H 2 -to- Bry ratio will be 
high, regardless of viewing angle. On the other hand, 
UV radiation from massive stars excites Bry emission as 
well as H 2 . Bry is often very bright, and when viewed 
face-on, one naturally finds relatively low H 2 -to- Bry ra¬ 
tios. In 30 Dor, the line ratio is no greater than 0.5 
(see § 13.41) , which favors fluorescence as the dominant H 2 


excitation mechanism. We note that in regions viewed 
edge-on where Bry is seen spatially offset from H 2 , such 
as the pillars in Area C, the H 2 -to- Bry ratio will be very 
low. Nevertheless, the presence of both Bry and H 2 ex¬ 
cludes shock excitation, regardless of viewing angle. 

5.2.1. H 2 and 8 / 1 m Emission Correlation 

The 8 /mi emission adds additional information for di¬ 
agnosing the origin of H 2 emission in 30 Dor. Emis¬ 
sion at 8 (im is largely dominated by Polycyclic Aro¬ 
matic Hydrocarbon (PAH) emission, which is excited by 
far-UV radiation in PDRs. Shocks, on the other hand, 
will d estroy the PAH molecules and suppress the emis¬ 
sion dFlower fe Pineau des Forets! f2003fc iMicelotta et al.1 
I2010T ). Therefore spatial and morphological correlations 
between the H 2 and 8 /im emission provide hints of 
the excitation mechanism of the observed H 2 emission. 
Good morphological correlations between the two distri¬ 
butions would suggest that H 2 emission is predominantly 
fluorescence excited, just like the 8 /mu emission; poor 
morphological correlations would otherwise imply that 
the FUV radiation is not a major excitation source of 
the observed H 2 emission. The only excessive H2 emis¬ 
sion related to shock activities is found near protostars 
at 4.5 /im, e.g. iCvganowski et all (|2008l) and lLee et al.1 

(IM1 . 

We convolved the H 2 image and Spitzer IRAC 8 /im 
image to the same spatial resolution of 2.0", and the su¬ 
perimposed image is shown in Figured The H 2 emission 
is displayed in red, while the 8 gm emission is in green. 
The overall morphologies of H 2 and 8 (im emission cor¬ 
relate very well in the entire 30 Dor nebula, as well as 
in all three areas of interest A, B, and C. In Figure E3 
we plotted the surface brightness of H 2 and 8 /iin emis¬ 
sion in Area A. The H 2 and 8 /im emission intensities 
also show a rather tight correlation. The good morpho¬ 
logical and intensity correlations between H 2 and 8 /im 
emission strongly implies that the observed H 2 emission 
is predominantly excited by the FUV radiation in the 
PDR. 

In several places within the nebula, especially Area C, 
H 2 and 8 /im emission highlight pillar structures. These 
are enclosed by thicker layers of ionized gas traced by 
Bry (noted in , which in turn have sharp boundaries 
that presumably arise from the pressure of the x-ray 
emitting hot gas. An example is shown in a zoom-in 
figure in Figure [H In other locations, such as the finger 
in Area A, the 8 /im emission shares an inner edge with 
the Bry emission, as though the grains responsible for 
this emission permeate the photoionized gas. We intend 
to further pursue the physical interpretation of these cor¬ 
relations in a future paper. 

5.3. Bright Emission Line Contamination 

We use our models to further explore the bright emis¬ 
sion line contamination issue in our H 2 image. Bry at 
2.17/rm and H 2 1-0 S(l) at 2.12 /im are present in the 
modeled continuum spectra and are very close to the cen¬ 
ter of the Ks band filter, as expected. In addition, there 
are He I 2 4 P-2 4 S emission line at 2.06 /im, H 2 1-0 S(2) 
at 2.03 gm, H 2 1-0 S(0) at 2.22 /im, and H 2 2-1 S(l) at 
2.25 /mu. Table [2] summarizes the emission lines in the 
Ks band. 
















































Table 2 

Emission lines within the Ks 
bandpass 
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Figure 5. Spatial and morphological correlations between H 2 and 
8 fim emission in 30 Doradus. H 2 emission is shown in red, and 
the image is convolved to a 2.0 // resolution to match that in the 
8 /im image. The 8 pm emission is shown in green. R136 is marked 
by a cyan cross. 




Figure 6. Surface brightness of H 2 and 8 pm in Area A of 30 
Doradus. The color bar indicates the number of pixels in each 
surface brightness bin of H 2 and 8 fi m data. The H 2 and 8 fi m 
intensities show an excellent correlation in this area, implying that 
the H 2 emission originated inside the PDRs. 


The H 2 1-0 S(0), H 2 1-0 S(2), and H 2 2 1 S(l) lines 
are negligible, for they contribute < 10% of the Bry flux 
in the Ks band. The He I line is the most plausible source 
of continuum contamination in our data. The calculated 
line ratios He I to Bry at densities < 10 4 cm -3 are 0.6 
to 0.7, which suggests that the He I line will contribute 
60% to 70% of the Bry intensity to the continuum level, 
in agreement with the empirically evaluated 61% level in 
§ l2.3.1l and the 70% level in M16 (LOO). Since Bry emis¬ 
sion is 7% of the Ks continuum, the He I contamination 
in the filter will be up to 4%, assuming the distribution 
is the same. 


6. SUMMARY 

We present the first and fully calibrated H 2 1-0 S(l) 
emission image of the entire 30 Doradus nebula, as well 
as a Bry image. In the data reduction process, we con¬ 
firmed Bry and He I emission line contamination in the 
Ks continuum via empirical analysis and Cloudy simula¬ 


Line 

Wavelength 

Br7 

2.17 

HeI 2 !p-2 !S 

2.06 fi m 

H 2 1-0 S(2) 

2.03 fim 

H 2 1-0 S(l) 

2.12 [im 

H 2 1-0 S(0) 

2.22 fim 

H 2 2-1 S(l) 

2.25 fin i 


tions. The error in contamination-corrected H 2 images 
is estimated to be ~ 4%. 

The overall morphology of H 2 correlates well with 
Bry and CO emission, implying that the observed H 2 
originates from the PDRs in 30 Doradus. The bright¬ 
est H 2 -emitting areas (Area A and Area B) are PDRs 
viewed face-on located behind the ionizing source R136, 
and the warm molecular clouds traced by H 2 appear to 
be clumpy. Those regions also trace the CO morphol¬ 
ogy well, further indicating that these PDRs are face-on, 
with layers of ionized gas, warm molecular gas, and cold 
molecular clouds. Discontinuity of H 2 and Bry mor¬ 
phology is found at the outer shells and pillar features 
(such as Area C) in 30 Dor , where the H 2 pillars are 
encompassed by Bry envelops of sharp boundaries. This 
suggests that we are viewing the the shells and pillars 
of the H II region edge-on. The mean projected separa¬ 
tion between the Bry envelope and H 2 clumps is 3 pc 
(Ay < 0.7), indicating that H 2 emission is formed in the 
PDRs close to the surface of the molecular clouds. The 
density of H 2 -emitting gas is inferred from the observed 
H 2 -to- Bry line ratios, the Cloudy model results, as well 
as theoretical predictions of H 2 emission in PDRs. The 
molecular gas density is estimated to be < 10 4 cm -3 . 
Low H 2 -to-Bry line ratios (< 0.5), as well as good mor¬ 
phological correlations between the H 2 and 8 /im emis¬ 
sion, implying again that the observed H 2 emission is 
excited by FUV radiation. 

While it requires multiple transitions of ro-vibrational 
H 2 lines to constrain the excitation mechanisms and 
physical parameters of the ISM, such data do not yet ex¬ 
ist. Our imaging observations suggest that the observed 
H 2 emission likely arises from a lower density layer of the 
PDR near IFs in 30 Doradus. We found no sign of shock- 
excited H 2 emission, and all indications were consistent 
with fluorescent excitation. 
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